Simulations of deposition growth models in various dimensions. 

Are overhangs important? 
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We present simulation results of deposition growth of surfaces in 2, 3 and 4 dimensions for ballistic 
deposition where overhangs are present, and for restricted solid on solid deposition where there are 
no overhangs. The values of the scaling exponents for the two models are found to be different, 
suggesting that they belong to different universality classes. 

PACS numbers: 68.55.Jk,05.70.Ln,61.50.Cj,Cd,64.60.Ht 



The deposition growth of surfaces JjJ has been a sub- 
ject of long continual theoretical and experimental inter- 
est [Q due to its relevance to non-equilibrium processes 
in general as well as its possible role in surface technol- 
ogy. The profile of the deposited surface gradually rough- 
ens under the stochastic accumulation of particles, and 
early simulations by Family et. al. Q suggested that 
the surface roughness exhibits a dynamical scaling be- 
haviour. That is, the height-height correlation function, 
G(r — r',t) = ([h(r,t) — h(r' ,t)] 2 )^ , scales with time, t, 
and separation, I — \r — r'\, as 



(1) 



h(r, t) is the height of the surface at position r and at time 
t. The dynamical scaling behaviour is characterised by 
the roughness exponent, a, and the dynamical exponent, 
(3, with z = a/ (3. The scaling function f(x) behaves as 
f(x) — x 13 for x <C 1 and f(x) — constant for x 3> 1. 
Thus, the surface roughness grows as G(t) ~ t 13 initially, 
independent of size, and for a given size, £, the roughness 
saturates after a sufficiently long time such that G(£) 
scales with £ only as G(£) ~ £ a . 

Numerous simulations in a variety of growth models 
(QilljlJ^l have since confirmed the hypothesis of dynam- 
ical scaling, including models which allow overhangs to 
form and models where overhangs are not allowed. An 
overhang is formed when a particle sticks at a position 
higher than the height of the surface at that point, such 
the space below the particle is not filled. Simulations of 
the restricted solid on solid model |l(]|ll| where incom- 
ing particles fall directly onto the surface such that no 
overhangs can form, and may only stick at a site if the re- 
sulting nearest neighbour height differences are less than 
some predetermined value, have led to a further consen- 
sus that the values of the scaling corresponds to that of 
the Kardar-Parisi-Zhang equation [[L2|, 
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at 



vV 2 h+^- (V/i) 2 +ry, 



was derived by assuming that the surface grows uniformly 
in the direction of the local normal. The exponents ob- 
tained are exact in 2 dimensions p2| , and numerically 
determined in higher dimensions MM- 

Results from simulations of the ballistic deposition 
model where incoming particles stick at the first 
point of contact and thus allow overhangs are more con- 
troversial. At present there is no clear consensus as to 
whether or not this system belongs to the same univer- 
sality class as that described by the Kardar-Parisi-Zhang 
equation |l4| , or whether the presence of overhangs leads 
to a different set of scaling exponents. Early results by 
Meakin et. al. gave a = 0.47 and (3 = 0.331 in 2 di- 
mensions, and a — 0.33 and (3 = 0.24 in 3 dimensions, in 
agreement with Kim and Kosterlitz's approximate for- 
mula |of« = 2/(d + 2) and (3 = l/{d + 1) for the 
Kardar-Parisi-Zhang equation. More recent results sug- 
gest that the values of the scaling exponents may, in fact, 
be different. Baiod et. al. obtained (3 = 1/3 in 2 di- 
mensions and a = 0.3 and (3 — 0.22 in 3 dimensions; 
off-lattice simulations have also given (3 — 0.343 in 2 
dimensions pi , but a clear scaling behaviour was not 
observed in 3 dimensions Jl7t . 

In this Letter, we report results of simulations of bal- 
listic deposition and restricted solid on solid growth. We 
find that the values of the scaling exponents for the bal- 
listic deposition model are different to those of the re- 
stricted solid on solid model. A summary of our results 
is given in Table 1. 

TABLE 1. Scaling exponents obtained from our simu- 
lations. 
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2 


0.45 0.32 


1.40 


0.50 0.33 1.50 


3 


0.26 0.21 


1.24 


0.40 0.25 1.60 


4 


- 0.12 




0.29 0.18 1.61 



where r\ is a random variable. This equation is believed 
to be a continuum description of deposition growth, and 
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introduce significant corrections. Instead, by rewriting 
equation ([!]) as 



FIG. 1. Ballistic deposition for 2 dimensions. The inset 
shows the plot of In G(£, t) versus log 2 I at the end of the 
simulation. 

We also find that while on-lattice simulations give ex- 
cellent scaling behaviour for the restricted solid on solid 
model, the same is not true for ballistic deposition. Quasi 
offdattice simulations were therefore carried out for the 
ballistic deposition model. Namely, each axis of a surface 
of size L d ~ Y particle diameters is divided into nL points 
such that incoming particles can be centred on any one of 
these points. For n large, the surface approaches a con- 
tinuum, and for n = 1, we recover the on-lattice model. 
The height of the surface at a position r is defined to 
be the height of a new particle if it fell onto the surface 
at r. We found that n — 3 is sufficient to give a good 
scaling behaviour, and no differences were found in the 
results with n = 5, 7 and 10. We performed simulations 
in 2, 3 and 4 dimensions for both models and the simula- 
tions are run until the equivalent of at least 2000 layers 
of atoms have been deposited. The total number of par- 
ticles deposited in each simulation is over 2 x 10 9 . The 
minimum time required for each run is twenty-four CPU 
hours on a DEC Alpha 400 workstation. To obtain good 
statistics, averages over many runs were often needed. 

In figure [j], the correlation function, G(£, t), for a 2 di- 
mensional ballistic deposition simulation is plotted versus 
time in a log-log plot. The largest system size considered 
is £ = 2 20 . For the larger values of £, the roughness has 
not saturated within the time scale of the simulation. In 
the dynamical scaling region, we see a clear power law 
behaviour, G(t) ~ t@. Also shown in the inset is a plot of 
lnG(£, t) versus log 2 t for the data at the end of our sim- 
ulation. For the smaller sizes where saturation has been 
reached, we also find a roughly linear dependence of In G 
on log 2 £ in agreement with the predictions of dynamical 
scaling. 

Direct extrapolation of the scaling exponents from the 
gradients in the log-log plots turned out to be difficult 
because cross-over effects due to the transition from the 
dynamical scaling regime to the saturated scaling regime 



\nG(£,t) - (3\nt = F(aln£ - /31ni), 
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we can obtained good estimates of the dynamical and 
roughness exponents by collapsing our data for all sizes 
and all times considered. This, in fact, provides a way of 
checking also whether the data corresponds to just one 
scaling regime, or whether there is also a cross-over be- 
tween different universality classes with different scaling 
exponents. We note that the surface roughness during 
the initial few time steps is strongly influenced by tran- 
sient effects, and have been discarded in the data col- 
lapse. 

The collapsed data for the 2 dimensional ballistic depo- 
sition result is shown in figure |^. Data for £ ranging from 
2 2 to 2 19 are used in the plot, with over 6 x 10 9 particles 
deposited. The values of the exponents used are a = 0.45 
and f3 = 0.32. We have also carried out simulations of 
the restricted solid on solid model in 2 dimensions, and 
found that a = 0.50 and j3 = 0.33, in agreement with the 
results of previous simulations. 



FIG. 2. Collapsed data for 2 dimensional ballistic deposi- 
tion simulation, with a = 0.45 and (3 = 0.32. The data used 
ranges from t = 2 2 to 2 19 . The inset show the collapse ob- 
tained with the exponents obtained for the restricted solid on 
solid model. 

The collapsed data for the 3 dimensional simulations 
are shown in figure | The upper diagram is for the 
restricted solid on solid model. The size of the system 
considered is 2 10 x 2 10 , and over 2 x 10 9 particles were 
deposited in a run. The data presented represents the 
average over seven independent runs, and include values 
for £ ranging from 2 2 to 2 9 . The values of the scaling 
exponents obtained in this case are a = 0.40 and (3 — 
0.25. This is in agreement with the approximate formula 
of Kim and Kosterlitz ||, but the value of (3 obtained 
is greater by 0.01 than that observed more recently by 
Ala-Nissila et. al. Ol. 
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Ballistic deposition simulations in 3 dimensions are 
also carried out for systems with size equal to 2 10 x 2 10 
particle diameters, with three subdivisions per particle 
diameter. The collapsed data are shown in the lower di- 
agram of figure H Again, over 2 x 10 9 particles were 
deposited per run, and the results presented represent 
the average over ten runs with data for I — 2 2 to 2 9 used 
in the data collapse. The simulations were also carried 
out with seven subdivisions per particle diameter and no 
difference was found. The value of the scaling exponents 
in this case are a = 0.26 and f3 = 0.21, significantly lower 
than the corresponding values for the restricted solid on 
solid model. 



the roughness exponent in 4 dimensions, however, give a 
value of a ~ 0.12. The uncertainty in this case is due 
to strong fluctuations in the roughness as a result of the 
small system size, and the number of runs required to 
obtain better statistics is prohibitively large. For the re- 
stricted solid on solid model, the fluctuations are smaller 
even in 4 dimensions and we have been able to obtain 
reliable values for both the dynamical and the roughness 
exponents. These are a = 0.29 and (3 = 0.18, in good 
agreement with those obtained by Ala-Nissila et. at [[ll| . 
Again, in accordance with the trend observed in lower di- 
mensions, the exponent of the ballistic deposition model 
is lower than that of the restricted solid on solid model. 



FIG. 3. Collapsed data for the 3 dimensional restricted 
solid on solid simulation (upper diagram), and the ballistic 
deposition simulation (lower diagram) for i = 2 2 to 2 9 . The 
inset in the lower diagram shows the collapse obtained if the 
exponents obtained from the restricted solid on solid model 
were used instead. 

We have also carried out simulations in higher dimen- 
sions for both models. However, due to computational 
difficulties, we are restricted to relatively small sizes. 
For the ballistic deposition model the largest size pos- 
sible in 4 dimensions or higher is still too small for the 
dynamical scaling regime to be observed. Estimates of 



FIG. 4. A cross-section of ballistic deposition growth in 3 
dimensions taken at a height of 200 particle diameters. The 
length of the horizontal and vertical axis correspond to 100 
particle diameters. 

We have found that variations in the value of either a 
or (3 by as little as 0.01 is sufficient to give clear deterio- 
ration of the data collapse plots. The values we present 
are therefore accurate to the figures quoted. The most 
important implication of this is that from our results, 
the dynamical scaling behaviour of the ballistic deposi- 
tion model and the restricted solid on solid model be- 
long to different universality classes. We have shown in 
the inset to figures || and || what happens when we try 
to collapse the ballistic deposition data with the expo- 
nents obtained from the corresponding restricted solid 
on solid simulations. It is clear from the diagrams that 
even in 2 dimensions where the differences between the 
values of the scaling exponents for the ballistic deposition 
model and those of the restricted solid on solid model are 
apparently small, a satisfactory data collapse cannot be 
obtained. In view of the belief that the dynamics of the 
restricted solid on solid model corresponds to that of the 
Kardar-Parisi-Zhang equation, our results would there- 
fore further suggest that the Kardar-Parisi-Zhang is not 
appropriate in describing deposition growth in situations 
where overhangs are dominant. Indeed, our values of the 
scaling exponents for the ballistic deposition model in 2, 
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3 and in 4 dimensions lie outside the range of the values 
for the Kardar-Parisi-Zhang exponents . 



ory for deposition growth in the presence of overhangs 
remain an important challenge. 
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FIG. 5. The fraction x of sites occupied along a substrate 
dimension versus the substrate dimension. 

We have also tried to examine the structure of the solid 
formed by growth under ballistic deposition conditions. 
Figure ^ show a cross-section of the bulk formed in a 3 
dimensional ballistic deposition simulations. The cross- 
section corresponds to a height of 200 particle diameters 
from the substrate, and is taken after all the particles at 
this height are covered. The cross section shown corre- 
sponds to an area of 100 x 100 particle diameters. We 
hnd that there are very few connected lines, and no con- 
nected rings in the cross section. In addition, we have 
calculated the fraction of sites, x, which are occupied in 
a linear direction from the average density, p. For a d s 
dimensional surface, the density is given by p = x ds . In 
figure U a plot of x versus substrate dimension is shown. 
The results indicate that of order 0.4 of the sites along a 
line on the surface are occupied in all dimensions. This, 
together with the cross section plot, corroborates with 
the idea that particles grow on the edges of overhangs, 
and almost immediately branch off to form a complex 
tree like structure. 

In summary, we have found that the presence of over- 
hangs is an important factor in determining the scaling 
properties of deposition growth. In a model such as bal- 
listic deposition, overhangs will form when the local sur- 
face gradient exceeds a critical value corresponding to the 
presence of a sharp step in the surface profile. In such a 
situation, the next particle will stick to increase the lat- 
eral size of the overhang region rather than to reduce the 
surface gradient by falling to the lower surface. Thus, as 
overhangs begin to form, they will tend to increase the 
lateral correlation at a fast rate, and the surface will no 
longer grow in the direction of its local gradient. The 
result may be an anisotropic growth which when coarse 
grained lead to broader and flatter structures. Although 
such a picture can give a behaviour consistent with the 
results of our simulations, the search for a proper the- 
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